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We compute the distribution of the work done in stretching a Gaussian polymer, made of N 
monomers, at a finite rate. For a one-dimensional polymer undergoing Rouse dynamics, the work 
distribution is a Gaussian and we explicitly compute the mean and width. The two cases where the 
polymer is stretched, either by constraining it's end or by constraining the force on it, are examined. 
We discuss connections to Jarzynski's equality and the fluctuation theorems. 
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I. INTRODUCTION. 

Classical thermodynamics does not give us all the de- 
tails about a nonequilibrium process. For example con- 
sider a nonequilibrium process during which we perform 
work W on a system kept in contact with a heat bath 
at some fixed temperature T . The system starts from an 
equilibrium state described by the temperature T and 
some other parameter, say A (e.g. volume). During the 
process the parameter changes from its initial value Xi to 
a final value A/. At the end of the process the system 
need not be in equilibrium but will eventually relax to 
an equilibrium state described by T and A / . The second 
law then tells us that 

W > AF (1) 

where AF is the difference of free energy between the 
two equilibrium states. The equality holds if the pro- 
cess is reversible. For an irreversible process what other 
information can one extract from a measurement of the 
work done ? First note that for specified initial and final 
values of the parameter A and a fixed path X(t) connect- 
ing them, the work done will not have a unique value. 
Every time we repeat the process we will get a different 
work done because: (a) the initial microscopic state we 
start from may be different and (b) for a given initial 
microscopic state the time evolution is not unique since 
the system is in contact with a heat bath. Thus we will 
get a probability distribution for the work done and it is 
of interest to examine the properties of this distribution. 
Recently there has been a lot of interest on issues related 
to properties of such distribution functions. Two very 
interesting results involving universal properties of these 
distributions have been proposed. 

The first is a very surprising exact equality obtained by 
Jarzynsk: 1 which states that: 

(e-W) = e~? AF (2) 

where the average is over the work distribution function. 
This result seems to hold very generally and should be 
compared with the inequality in Eq. ^ that one obtains 
from usual thermodynamics. 

The second set of results are obtained when one looks 
at the probability distribution of various non-equilibrium 



quantities (including W) such as, for example, the en- 
tropy production. In this case some new fluctuation theo- 
rems have been proposed 2 3 4 . These theorems were origi- 
nally derived for deterministic systems but have also been 
proved for stochastic systems 5 ' 6 ' 7 " These theorems look at 
the ratio of the probabilities of positive to negative en- 
tropy production during a nonequilibrium process and 
thus give some measure of "second law violations" which 
can be significant if one is looking either at small systems 
or at small time intervals. There are two versions of the 
fluctuation theorem, the steady state fluctuation theorem 
(SSFT) and the transient fluctuation theorem (TFT). In 
the former case one looks at a system in a nonequilib- 
rium steady state and the average entropy production 
rate is examined. In the transient fluctuation theorem a 
system is initially prepared in thermal equilibrium and 
one looks at the entropy produced in a finite time r. 
An important point to note is that the definition of en- 
tropy production in small (non-thermodynamic) systems 
and in a nonequilibrium situation is somewhat ad-hoc 
and various definitions have been used. A number of 
authora&iLi2iiiii2ii2have looked, both theoretically and in 
experiments, at fluctuations of quantities such as work, 
power flux, heat absorbed, etc. during a nonequilibrium 
process. In an interesting work it was shown by Crooks 6 
that the Jarzynski equality and the TFT are connected. 

Finding universal properties of various nonequilibrium 
distribution functions is of obvious interest. At the same 
time the explicit forms of the distributions for different 
systems is clearly of interest too. Infact for systems such 
as polymers these are experimentally accessible^ and it 
seems plausible that they can give information on the 
dynamics of the system. There has not been much work 
in this direction. In a recent paper— Speck and Seifert 
have shown that in the limit of slow driving the work 
distribution becomes a Gaussian. They consider systems 
(which could be nonlinear) evolving through Langevin 
dynamics. Apart from this work, most other explicit cal- 
culations of nonequilibrium distribution functions have 
considered single particle systems. 

In this paper we consider the well-known model of a 
flexible polymer whose motion is governed by Rouse dy- 
namics. We look at the work done when the polymer is 
stretched at a finite rate. The work distribution func- 
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tions in different ensembles ( constant force and constant 
extension) are computed explicitly and the dependence 
of the distributions on switching rates and system sizes 
is examined. We discuss our results in the context of 
Jarzynski's equality and the fluctuation theorems (TFT). 



II. DISTRIBUTION OF WORK IN THE 
CONSTANT EXTENSION ENSEMBLE. 

We consider a one-dimensional Gaussian polymer 
whose energy is given by: 



(3) 



Z=1,ZV+1 



with j/o = and Un+i = c*(t) which is a specified function 
of time. For the moment we take the spring constants ki 
to be arbitrary. We assume the following Rouse dynamics 
for the chain: 

m = --(yi-yi-i)+—(yi+i-yi)+m l = l,2...N (4) 

7 7 



Since W is linear in y(0) and r\ both of which are Gaussian 
variables, it follows that the distribution of W will also 
be Gaussian. We then only need to find the mean and 
the second moment which we now obtain. We first state 
a few results on equilibrium properties of the Gaussian 
chain. The equilibrium free energy of the chain is given 
by 

Z(a) = J dyidy2---dyisfe~^ H where 
N+l , 



B = £?(w-w-i) 2 



1=1 

^y T Ay - h T y + ^k N+1 a 2 



(8) 



This leads to the following equilibrium free-energy of the 
polymer (apart from a-independent constants). 

F(a) = hn{Z) = h N+1 a 2 - Va^ (9) 
= \ka 2 



which can be written in matrix notation as 

dy 1 . 1 , , . 
JL = —Ay + ~h{t) + r ] 
at 7 7 



(5) 



where y T = {y x , y 2 — -Vn}, if = {vuVi, —Vn}, h T = 
{0,0, ...kN+itt{t)} and the noise satisfies (i]i(t)) = and 
(m{t)Vm{t')) = 2/(pj)8 lm 8(t - t') . The matrix A is 
tridiagonal with elements A;^ = (fc; + fcz+i), A-i.i+i — 
—ki+i, A; j_i = — h\. The general solution of this equa- 
tion is given by: 



y(t) = G(t)y(0)+ [ dt'Git-t 
Jo 



,Mt>) 
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dt'G(t-t')r}(t') 
(6) 

where G(t) = e At /i. The work done, as defined by 
Jarzynski, is then given by: 

Wj = f ^adt 
Jo oat 



= k N+ i I (a- y N )adt 
= ^(a 2 (r)-a 2 (0))^ / ,////'■,,. 



This work is equal to the true mechanical work done by 
the external force which we will call W (thus in this case 
W = Wj). We now plug in the solution for y^ from 
Eq.EJto get 

W = ^V(r) -a 2 (0)) - f T dthT[G(t)y(0) 



f dt'G(t - 1') 
Jo 
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dt'G(t - t')r)(t')\ (7) 



where 1/fc = 1/ki + 1/&2 + •••1/fejv+i. The mean posi- 
tions of the particles and their fluctuations can be easily 
computed at any force and are given by: 



(y) = a- 1 /* 



(10) 



(((y-(y))(y T ~(y T ))) = ^A- 1 (11) 

Mean and fluctuations of the work done : 

From Eq. 0we get for the mean work done: 

(W) = ^V(r) - a 2 (0)) £ dth T [G(t)(y(0)) 



- [ dt [ dt'h T G(t - t')h(t') 
7 Jo Jo 



(12) 



We do integration by parts so as to express everything 
in terms of the rates h. Using the equilibrium results in 
Eq. 19I10|I we finally get: 

(W)= F(a(r)) - F(a(0)) 

+ [ dt [ dt'h T (t)A- 1 G(t-t')h(l/)(13) 
Jo Jo 

The fluctuations of the work {(W — (W)) 2 ) = a 2 is given 
by: 



a 2 = 



T pT 



dt / dt'k T (t)G(t)([y(o) - <y(o))Mo) - ( y (o))] T ) 



xG(t')h(t')+ dti dt\ dt 2 di/ 2 h T (ti)G(ti - 1[) 
Jo Jo Jo Jo 

x(r 7 (tih T (4))G(t 2 -4)/ l (t 2 ) 




FIG. 1: Distributions of the work done in pulling a short 
polymer (N = 1) at different rates r = tr/t. 
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FIG. 2: Distributions of the work done in pulling a long poly- 
mer (iV = 100) at different rates r — tr/t. 



Using Eq. (fTT|) and the relation (rj(t)r] T (t')) = 
2/(/?7)£(i - t')\ this simplifies to: 



- [ dt [ dl/h T (t)G(t)A- 1 G(t')h(t') 
P Jo Jo 
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dti I dt 2 dt' 2 h T {t^G^ -t' 2 )G(t 2 -t' 2 )h(t 2 



Finally using the relation J* 2 dt' 2 G(ti — t' 2 )G(t2 — t 2 ) = 
^[G(h - t a ) - G(ti)G(ta)] we get 



<t 2 = - dt dt'h T (t)A- x G(t-t')h(t') (14) 
P Jo Jo 

The distribution of work done is then: 



P{W) = 



1 



(15) 



(27T ( 7 2 ) 1 /2 

As expected for a Gaussian process we find that 

(W) - AF = (3a 2 /2. (16) 

Note that in the present polymer model both the equi- 
librium free energy and the average work are indepen- 
dent of temperature while the width of the distribution 
a 2 depends linearly on temperature. It is easily ver- 
ified that the work distribution satisfies the Jarzynski 



equality (e 



-pw\ __ 



£ AF . On the other hand, with 



the present definition of work, the fluctuation theorem 
is not satisfied. However if we define the "dissipated 
work" Wdiss = W — AF then the distribution of Wdiss, 
P{Wdiss) satisfies the fluctuation theorem: 



P{W, 



cits 



P(-W dlss ) 



pw dis 



Since Wdiss gives some measure of deviation from a qua- 
sistatic and adiabatic process it seems reasonable to think 
of it as an entropy production term which is what usually 
appears in the fluctuation theorems. 



Some special cases: Let us consider the case where 
all spring constants are equal ki = k and the let us assume 
that the polymer is pulled at a constant rate so that 
a(t) = at/r. The effective spring constant is k/(N + 1) 

From 



and so the free energy of the polymer is 
Eq. ^]we get the spread in the work done: 



2(W+1)- 



2a 2 7 



[k 2 A' 



k 



7 fc 3 A-'V 



Ar/7 



1)}nn (17) 



where [...]jviv denotes a matrix element. The average 
work can be obtained from Eq. ^j] In the two limits of 
very slow (r — > oo) and very fast (r — > 0) processes, we 
get: 



2fc 2 a 2 7A^ 2 JV 1 Aa 2 jN 



ka 2 N 

^vTT) 



/3tt 2 t 
(Fast) 



(Slow) 



(18) 
(19) 



For an instantaneous pulling process the work done is 
simply W = fcjv+i [(a(r) - y N ) 2 - (a(0) - Vn) 2 ]/^ and 
the result in Eq. EH can t> e directly obtained. 

It is instructive to plot the work distributions for differ- 
ent pulling rates. To see the effect of the polymer length 
on the distribution we consider two cases: (i) A short 
polymer with N = 1, k = 1, 7 = 0.1, a = 1, (3=1 and 
(ii) a long polymer with N = 100, k = 1, 7 = 0.1, a = 
10, (3=1 (in arbitrary units). In each case the param- 
eter values are chosen so that the change in equilibrium 
free energy given by AF = 2{ t ^ +l \ ) is of the order of 1/(3 . 
The pulling rate has to be compared with the relaxation 
time of the polymer which is given by tr = 7/ X sm where 
X sm = Ak sin ( 2(N+i) ) ^ s ^ ne smallest eigenvalue of A. 
For large N we get tr = 7 TV 2 / (kir 2 ) . For the cases (i) 
and (ii) we numerically evaluate Eq. 1171 for pulling rates 
r = tr/t = 0.1, 1, 10, 100. The resulting distributions 
are plotted in Fig.^ an d Fig- 121 F° r the long polymer the 
probability of negative work realizations is quite small. 
We can increase their probability by increasing the tern- 
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perature which broadens the distributions while keeping 
the mean unchanged. 



III. DISTRIBUTION OF WORK IN THE 
CONSTANT-FORCE ENSEMBLE. 

Next we compute the work distribution in the 
constant-force ensemble. Instead of constraining the end 
of the polymer we apply a time-dependent force f(t) on 
it. The time-dependent Hamiltonian of the system is now 
given by: 



H 



(20) 



1=1, N 



and the equations of motion are again: 

dy 1 A 1 , , v 
-f = —Ay + -h(t)+ V 
at j j 

with h T (t) — [0, 0, ...f(t)]. In this case we note that 
the the generalized work Wj, as defined by Jarzynski, 
is given by Wj = — J Q r dty^f{t)dt and is not equal to 
the true mechanical work done on the system which is 
W = Jq dtf{t)y N . In this paper we will compute the 
distribution of the true mechanical work W. We again 
find that W has a Gaussian distribution with the follow- 
ing mean and variance: 

(W) = - ( dth T (t)h(t) - i f dt [ dt'h T (t) 
1 Jo 1 Jo Jo 

i r 

xG(t - t')Ah(t') - - / dth T (t)G(t)h(0) 
9 Jo 

a 2 = A f dth T (t)h(t) - ^ / dt f dt'h T (t) 
P7 Jo PI Jo Jo 

xG(t - t')Ah(t') 

For /(0) = 0, we get (W) = fcr 2 which means that P(W) 
satisfies the fluctuation theorem. It is then natural to 
again ask if W is some measure of entropy production. 
If this was so then W should vanish for an adiabatic 
process. To check this we first express (W) in terms of 
the rate hit). We get 

(W) = ^(tJA-^Ct) - /i T (0)A _1 /i(0)] 

-h T {r) [ dtA- 1 G(T-t)h(t)+ [ dt [ dt'h T (t) 
Jo Jo Jo 

xA- l G{t-t')h{t') 



Hence for an adiabatic process we get (W) — — AG(/) 
where G(/) = — h T A~ 1 h/2 is the polymer free energy in 
the constant force ensemble. Thus in this case (W) is not 
zero for an adiabatic process and so it is not an obvious 
measure of entropy production even though it does satisfy 
the fluctuation theorem. 

Interestingly since a 2 = 2(W)//3, thus even for an 
adiabatic process, the work-distribution does not tend to 
a S function as one might naively expect. However in the 
thermodynamic limit N — > oo the width approaches zero 
as a ~ 1/N 1 / 2 so in this limit the usual expectation is 
indeed satisfied. 



IV. CONCLUSIONS. 

In conclusion, in this paper we have computed explic- 
itly the distribution of work done when a polymer is 
stretched at a finite rate. We examine different ensem- 
bles and look at different definitions of work. As has been 
noted by earlier authors, for different definitions of the 
work, the corresponding distributions can have quite dif- 
ferent properties. In the constant extension ensemble the 
generalized work Wj is the same as the true mechanical 
work W and Jarzynski's identity is satisfied. In this case, 
the fluctuation theorem is satisfied by a different quantity 
Wdiss which does seem like a quantity which gives some 
measure of entropy production. In the constant force en- 
semble, Wj is different from the true work W which does 
not satisfy the Jarzynski identity. On the other hand the 
distribution of W satisfies a fluctuation-theorem like re- 
lation. However in this case we find that it is not possible 
to identify the work as a measure of entropy production. 

As a practical use, the Jarzynski identity has been pro- 
posed as an efficient method for computing equilibrium 
free energy profiles from noncquilibrium measurements, 
both in simulations and experimentsASiii. For the speci- 
cific case of polymers, we expect that a combination of 
simulations and our exact results on the work distribu- 
tion, should lead to better estimates on efficiency and of 
errorsi&ASiSI (and how they depend on rates and system 
sizes) involved while using the nonequilibrium methods 
in free energy computations. 

For non-Gaussian models of polymers such as, for ex- 
ample, semiflexible polymers, the work distribution func- 
tion is likely to be non-Gaussian. For small pulling rates 
one again expects a Gaussian distribution. It will be 
interesting to compute such distributions explicitly and 
study their dependences on rates, system sizes and other 
parameters such as the rigidity of the polymer. 
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